Brane universes tested by supernovae 



Mariusz P. DabrowskiQ 

Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland 

Wlodzimierz Godlowskijj] and Marek Szydlowski[] 

Astronomical Observatory, Jagiellonian University, 30-244 Krakow, ul. Orla 171, Poland 

(Dated: February 5, 2008) 

We discuss observational constrains coming from supernovae la || imposed on the behaviour of 
the Randall-Sundrum models. In the case of dust matter on the brane, the difference between the 
best-fit general relativistic model with a A-term j| and the best-fit brane models becomes detectable 
for redshifts z > 0.6. It is interesting that brane models predict brightergalaxies for such redshifts 
which is in agreement with the measurement of the z = 1.7 supernova JfJ and with the New Data 
from the High Z Supernovae Search Team We also demonstrate that the fit to supernovae data 
can also be obtained, if we admit the "super-negative" dark energy p — — (4/3)g on the brane, where 
the dark energy in a way mimics the influence of the cosmological constant. It also appears that the 
dark energy enlarges the age of the universe which is demanded in cosmology. Finally, we propose 
to check for dark radiation and brane tension by the application of the angular diameter of galaxies 
minimum value test. 

PACS numbers: 04.20. Jb,04.65.+e,98.80.Hw 



The idea of brane universes has originated from Hof ava 
and Witten followed by Randall and Sundrum ||. 
Brane models admit new parameters which are not 
present in standard cosmology (brane tension A and dark 
radiation IA). From the astronomical observations of su- 
pernovae la |3|, U one knows that the universe is now ac- 
celerating and the best-fit model is for the 4-dimensional 
cosmological constant density parameter S7a (4) ,o = 0.72 
and for the dust density parameter f2 OT ,o = 0.28 (in- 
dex "0" refers to the present moment of time). In other 
words, only the exotic (negative pressure) matter in stan- 
dard cosmology can lead to this global effect. On the 
other hand, in brane models the g 2 quadratic contribu- 
tion in the energy density even for a small negative pres- 
sure, contributes effectively as the positive pressure, and 
makes brane models less accelerating. In this paper we 
argue that in order to avoid this problem one requires 
much stronger negative pressure p < — g matter to be 
present on the brane (cf. Ref. ||) unless the new HZT 
data shows the z > 1 supernovae are brighter than ex- 
pected H 7|. 

In Ref. [E ] we gave the formalism to express dynamical 
equations in terms of dimensionless observational density 
parameters fl. Following Refs. || [l0| |T^] we introduce 
the notation useful for this purpose. In this notation the 
Friedmann equation for brane universes takes the form 
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where a(t) is the scale factor, k = 0,±1 the curvature 
index, Au\ the 4-dimensional cosmological contant, and 
7 the barotropic index (p = (7 — 1)^, p - the pressure, g 
the energy density) . In Eq. ([!]) we have defined the ap- 
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is a of general relativistic nature, C\ comes as contribu- 
tion from brane tension A, and Cu as a contribution from 
dark radiation. Though in Refs. || [l2j the Eq. ([[]) was 
studied using qualitative methods we just mention here 
that the cases 7 = (cosmological constant), 7=1/3 
(domain walls) and 7 = 2/3 (cosmic strings) can be ex- 
actly integrable in terms of elementary |l(| or elliptic jj| 
functions. For other values of 7 = 4/3; 1; 2, the terms 
of the type l/o 8 and 1/a 12 appear, and the integration 
involves hyperelliptic functions. In particular, oscillating 
non-singular solutions appear for dark energy 7 = —1/3 
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In order to study observational tests we now define 
dimensionless observational density parameters jl0| 
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where the Hubble parameter H — a/ a, and the decel- 
eration parameter q = —da/ a 2 , so that the Friedmann 
equation ([!]) can be written down in the form 

n G R + + n k + n A(4) + n u = 1. (3) 

Note that flu m despite standard radiation term, can 
either be positive or negative. Using (||), the equation (|l|) 
can now be rewritten as (compare Eq.(10) of fLOf) 
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It is also useful to express the curvature of spatial sections 
by observational parameters by using (||) and (^) 

- fl k = y Qgr + + 2Skt-q-l- (5) 
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These relations (|j) and (||) allow to write down an ex- 
plicit redshift-magnitude relation (a generalized Hubble 
law) for the brane models to study their compatibility 
with astronomical data which is the subject of the present 
paper. Obviously, the luminosity of galaxies depends on 
the present densities of the different components of mat- 
ter content Q given by (0) and their equations of state 
reflected by the value of the barotropic index 7. 

Let us consider an observer located at r = at the 
moment t = to which receives a light ray emitted at t = t\ 
from the source of the absolute luminosity L located at 
the radial distance r\. The redshift z of the source is 
related to the scale factor a(t) at the two moments of 
evolution by 1 + z = a(to) / a(ti) = a Q /a. If the apparent 
luminosity of the source as measured by the observer is I, 
then the luminosity distance dh of the source is defined 
by the relation 



I = 
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and T>l is the dimensionless luminosity distance. The 
observed and absolute luminosities are defined in terms 
of K-corrected apparent and absolute magnitudes m and 
M. When written in terms of m and M, Eq.(pf) yields 



m(z) = M + 5log w [V L (z)}, 
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where M. = M — 51og 10 Hq + 25. For homogeneous and 
isotropic Friedmann models one gets [O 
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where S(x) = sinx for /C = — f2fe,o; S(\) — X for /C = 1; 
S(x) = smn f° r ^ = ^fe,0- From the Friedmann equation 
([j]) and the form of the FRW metric we have 
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Firstly, we will study the case 7=1 (dust on the brane; 
we will label SIgr by fl m )- The case 7 = 2/3 (cosmic 
strings on the brane) has recently been studied in Q 
where, in fact, flu and fl\ were neglected and where the 
term r2 mj o(l + z ) 3 was introduced in order to admit dust 
matter on the brane. This case was already presented in a 
different framework in Ref. |llj] . Secondly, we will study 
the case 7 = — 1/3 (dark energy on the brane Q - we will 
label this type of matter with fid instead of flan). 



Now we test brane models using the sample of Ref. [|| . 
In order to avoid any possible selection effects, we use 
the full sample (usually, one excludes two data points as 
outliers and another two points, presumably reddened, 
from the full sample of 60 supcrnovae). It moans that 
our basic sample is the sample A of Ref. || . We test our 
model using the likelihood method 

First of all, we estimated the value of M. from the 
sample of 18 low redshift supernovae, also testing our 
result by the full sample of 60 supernovae taking £l\ = 
0. We obtained M. = —3.39 which is in a very good 
agreement with the results of Refs. |ll| 14 1. Also, we 
obtained for the model of Ref. B same value of % 2 = 
96.5. 

Neglecting dark radiation f^o = we formally got 
the best fit (x 2 = 94.6) for Q k<0 = -0.9, Q mfi = 0.59, 
^a,o = 0.04, CIa,o — 1-27, which is completely unrealistic, 
because fl m ,o = 0.59 is too large in comparison with the 
observational limit (also fijt,o = —0.9 is not very realistic 
from the observational point of view). 

However, we should note that, in fact, we have an ellip- 
soid of admissible models in a 3D parameter space f2 TO ,o> 
^a,0i ^A( 4 ),o at hand. Then, we have more freedom than 
in the case of analysis of Ref. || where they had only an 
ellipse in a 2D parameter space flmfi, f2A (4) .o- For a flat 
model flkfi = we obtain "corridors" of possible mod- 
els (0). Formally, the best-fit flat model is fl m .o — 0.01, 
Q\.o = 0.09 x 2 = 94.7 which is again unrealistic. In the 
realistic case we obtain for a flat model £l m ,o = 0.25, 
n Xfi = 0.02, n A(4) ,o = 0.73 with x 2 = 95.6. One should 
note that all realistic brane models require also the pres- 
ence of the positive 4-dimensional cosmological constant 

(n Aw ,o~o.7). 

There is a question if we could fit a model with neg- 
ative £l\fi7 For instance, in a flat Universe we could fit 
the model with fl m .o = 0.35 (too much in comparison 
with the observational limits on the mass of the clus- 
ter of galaxies) f2 A ,o = -0.01, f2 A ,o = 0.66 (x 2 = 96.3). 
However, it is not possible to fit any models with more 
negative values of fl\,o, regardless the Universe is flat or 
not. 

In Figjl] we present plots of redshift-magnitude rela- 
tions against the supernovae data. One can observe that 
in both cases (best-fit and best-fit flat models) the differ- 
ence between brane models and a pure flat model with 



is largest for 0.6 < z < 0.7 while it sig- 



nificantly decreases for the higher redshifts. It gives us 
a possibility to discriminate between general relativis- 
tic models and brane models when the data from high- 
redshift supernovae z > 1 is available. On the other hand, 
the difference between the best-fit general relativistic 
model with a A-term Q and the best-fit brane models 
becomes detectable for redshifts z > 0.6. It is interest- 
ing that brane models predict brighter galaxies for such 
redshifts which is in agreement with the measurement of 
the z = 1.7 supernova fa] and with the New Data from 
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FIG. 1: Redshift-magnitude relations for 7 = 1 brane uni- 
verses (dust on the brane). The top line is the best-fit flat 
model of Ref. || with fi m ,o = 0.28, fi A(4) ,o = 0.72. The bot- 
tom line is a pure flat model with fiA (4) ,o = 0. Between these 
two there are brane models with Q\ : o 7^ 0: lower - the best-fit 
model; higher - the best-fit flat model. 



the High Z Supernovae Search Team |7]]. In other words, 
if the farthest z > 1 supernovae were brighter, the brane 
universes would be the reality. 

One should note that we made our analysis without ex- 
cluding any supernovae from the sample. However, from 
the formal point of view, when we analyze the full sample 
A, all models should be rejected even on the confidence 
level of 0.99. One of the reason could be the fact that the 
assumed error bars are too small. However, in majority 
of papers another solution is suggested. Usually, one ex- 
cludes 2 supernovae as outliers, and 2 as reddened from 
the sample of 42 high-redshift supernovae and eventu- 
ally 2 outliers from the sample of 18 low-redshift super- 
novae. We decided to use the full sample A as our basic 
sample because a rejection of any supernovae from the 
sample can be the source of lack of control for selection 
effects. However, for completness, we also made our anal- 
ysis using samples B and C. It emerged that it does not 
significantly changes our results, though, increases qual- 
ity of the fit. Formally, the best fit for the sample B is 
(56 supernovae) (x 2 = 57.3): flk,o = —0.1 fl m ,o = 0.17, 
^a.o = 0.06, Oa (4)i o = 0.87. For the flat model we obtain 
(X 2 = 57.3): fl mfi = 0.12, fl Xfi = 0.06, fi A(1) , = 0.82, 
while for a "realistic" model (O m ,o = 0.25, fl\,o = 0.02) 
X 2 = 57.6. Formally, the best fit for the sample C (54 su- 
pernovae) (x 2 = 53.5) gives flk,o = (flat) fl m ,o = 0.21, 
^a,o = 0.04, fiA (4) ,o = 0-75, while for "realistic" model 
{fl mfi = 0.27, fl x ,o = 0.02) x 2 = 53.6. 

One should note that we have also separately estimated 
the value of M. for the sample B and C. We obtained 
A4 = —3.42 which is again in a good agreement with 
the results of Ref. |L3| (for a "combined" sample one 
obtains M. — —3.45). However, if we use this value in our 
analysis it does not change significantly the results (x 2 
does not change more than 1 which is a marginal effect 



for x 2 distribution for 53 or 55 degrees of freedom) . 
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FIG. 2: A redshift-magnitude relation for 7 = —1/3 brane 
universes (dark energy on the brane). The top and bottom 
lines are same as in Fig. |l|. The brane dark energy models 
plots are very close to the top line of Ref. || . 

In Fig. U we present a redshift-magnitude relation ( |l0| ) 
for brane models with dark energy (7 = —1/3). Note 
that the theoretical curves are very close to that of 
which means that the dark energy cancels the positive- 
pressure influence of the g 2 term and can simulate the 
negative-pressure influence of the cosmological constant 
to cause cosmic acceleration. From the formal point of 
view the best fit is (x 2 — 95.4) for flf.fi — 0.2, Q^o = 
0.7, fl\,o = —0.1, flu = 0.2, ^a (4) ,o = which means 
that the cosmological constant must necessarily vanish. 
From this result we can conclude that the dark energy 
p = —(4/3)f> can mimic the contribution from the A( 4 )- 
term in standard models. 

For the best-fit flat model (flk,a = 0) we have (x 2 — 
95.4): fl dfi = 0.2, fl K0 = -0.1, fl u = 0.2, fl K „ = 0.7. 
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FIG. 3: The age of the universe to in units of Hq 1 for the 
brane models with dust (0 < fl m ,o < 1 on the horizontal 
axis). Here flu,o = flk,o = 0, S1a,o = 0,0.05,0.1 (top, middle, 
bottom). 

Now let us briefly discuss the effect of brane parame- 
ters and dark energy onto the age of the universe which 
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according to (|l|) is given by 

H t = f {n GRfi x-^ +i + n x , oX -^ +i 

Jo 
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+fl u ,o + Qk,ox 2 + ^A (4) ,o^ 4 } 2 xdx, (11) 



where x — a/ag. We made a plot for the dust 7=1 
on the brane in Fig.|| which shows that the effect of 
quadratic in energy density term represented by fix is to 
lower significantly the age of the universe. The problem 
can be avoided, if we accept the dark energy 7 = — 1/3 || 
on the brane, since the dark energy has a very strong in- 
fluence to increase the age (jl6| . In Fig. [| we made a plot 
for this case which shows how the dark energy enlarges 
the age. 
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FIG. 4: The age of the universe to in units of Hq for the 
brane models with dark energy on the brane (0 < fld,o < 1 
on the horizontal axis). Here Qu,o = 0.2, fl\.o = 0.05,0 (top, 
bottom) which shows weaker influence of the brane effects to 
increase the age. 



Finally, let us study the angular diameter test for brane 
universes. The angular diameter of a galaxy is defined 
by 



d{z + lf 
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where d is a linear size of the galaxy. In a flat dust (7 = 1) 
univer se 9 has the minimum value z m j n = 5/4. It is par- 
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The angular diameter 6 for fix — 0.1, fir, 
= 0.72, and the two values of flu = 0.1, 
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(top, middle) in comparison with the model of Ref. j| with 



fin 



0.28, n A(4) =0.72 (bottom). 



value of flu = —1/3 (z„ 



2) in order to observe the 



minimum. 



More detailed analytic and numerical studies |T(| show 
that the increase of z m i n is even more sensitive to nega- 
tive values of fix (f? 2 contribution). Similarly as for the 
dark radiation fly, the minimum disappears for some 
large negative fix- Positive flu and fix make z min de- 
crease. In Fig. H we present a plot from which one can 
see the sensitivity of z m i n to flu- We have also checked 
|l6| that the dark energy fid has very little influence onto 
the value of z m in- 
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ticu larly interesting to notice that for flat brane models 



with fl\ « 0, ^A (4) ~ the dark radiation can enlarge the 
minimum value of 8 while the ordinary radiation lowers 
this value |lO|L i.e., 
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(13) 



for flu < 0. This is a general influence of negative dark 
radiation onto the angular diameter size for brane mod- 
els. One can also notice that there exists a restriction on 
the amount of negative dark radiation coming from ( p^ ) 
(flu > —1/3) which can serve as a test for the admissible 
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